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LOCALIZATION ERRORS IN SOLVING STOCHASTIC PARTIAL 
DIFFERENTIAL EQUATIONS IN THE WHOLE SPACE 

MATE GERENCSER AND ISTVAN GYONGY 


Abstract. Cauchy problems with SPDEs on the whole space are localized to 
Cauchy problems on a ball of radius R. This localization reduces various kinds 
of spatial approximation schemes to finite dimensional problems. The error is 
shown to be exponentially small. As an application, a numerical scheme is pre¬ 
sented which combines the localization and the space and time discretization, 
and thus is fully implementable. 


1. Introduction 

Parabolic, possibly degenerate, linear stochastic partial differential equations 
(SPDEs) are considered. In applications such equations are often given on the whole 
space, which makes the direct implementation of discretization methods problem¬ 
atic. Finite element methods, see e.g., m, m. na, m. [25] and their references, 
mostly treat problems on bounded domains and often under strong restrictions 
on the differential operators, denoted by L and M below. For finite difference 
schemes convergence results are available on the whole space, see e.g. [26] for the 
non-degenerate and m. m for the degenerate case, but the schemes themselves are 
infinite systems of equations. A natural way to overcome this difficulty is to “cut 
off” the equation outside of a large ball of radius R. A similar approach to obtain 
error estimates for a truncated terminal condition in deterministic PDFs of optimal 
stopping problems is used in [23]. 

The main results of the present paper are Theorems 12.61 and 13.61 Theorem 
12.61 compares solutions of two SPDEs whose data agrees on a ball of radius R and 
establishes an error estimate of order in the supremum norm. The proof relies 

on transforming fully degenerate SPDEs to zero order equations via the method of 
characteristics, whose analysis goes back to m and [18], see also the recent work 
[21] and the references therein. Once such a result is used to estimate the difference 
between the original equation and its truncation, one can approximate the truncated 
equation with known numerical schemes. Our choice is the finite difference method 
for the spatial and the implicit Euler method for the temporal approximation. The 
analysis of the former is invoked from [7] , while for the latter one requires an error 
estimate for the time discretization of the finite difference scheme, which is at this 
point a finite dimensional SDE. Of course this estimate needs to be independent of 
the spatial mesh size, and this is established in Theorem l3.101 To the authors’ best 
knowledge such an analysis of a full discretization is also a new result for degenerate 
SPDEs, and in fact even in the deterministic case it improves the results of [2] in 
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that we are not restricted to finite difference schemes which are monotone. An 
error estimate for the approximations obtained by localized and fully discretized 
SPDEs is given in Theorem 13.51 which is a special case of Theorem 13.61 where the 
accuracy of the approximations is improved by Richardson extrapolation in the 
spatial discretization. 

We mention that an alternative method for localization is to introduce artificial 
Dirichlet boundary conditions on a large ball, this approach is used for deterministic 
PDEs in, for example, in m, m, and to some extent, in [53]. The order of 
error in these works is e~^^, and so the method of the present paper yields an 
improvement even in the deterministic case. In the present context the method 
of artificial boundary condition would present additional issues. For instance, if 
the original SPDE is degenerate, then after introducing the boundary conditions 
the resulting equation may not even have a solution. Even if we do assume non¬ 
degeneracy, in the generality considered here - which is justified and motivated 
by the filtering equation in signal-observation models - there is no approximation 
theory of SPDEs on bounded domains with Dirichlet boundary condition, and 
therefore such a localization method does not help in finding an efficient numerical 
scheme. One indication of the problem is that regardless of the smoothness of the 
data and of the boundary, solutions of Dirichlet problems for SPDEs in general do 
not have continuous second derivatives, see HZ). 

Let us introduce some notation used throughout the paper. All random ele¬ 
ments will be given on a fixed probability space (fl, J-, P), equipped with a filtration 
{J't)t>o of cr-fields Pt C P. We suppose that this probability space carries a se¬ 
quence of independent Wiener processes adapted to the filtration {Pt)t>o, 

such that Wf — is independent of Pg for each k and any 0 < s < t. It is assumed 
that pQ contains all P-null subsets of If, so that (O, P, P) is a complete probability 
space and the tr-fields Pt are complete. By V we denote the predictable tr-field 
of subsets of D X [0,oo) generated by (J^t)t>o. We use the shorthand notation 
E“A = [E(A)]“. 

For p € [2,oo) and t? € M we denote by ,T-L) the space of measurable 

mappings / from into a separable Hilbert space 'H, such that 

|/| lp ,„ = ( / 1(1 + dxf'^ < oo. 

./R'i 

We do not include the symbol Ti, in the notation of the norm in Ti). Which 

% is involved will be clear from the context. We do the same in other similar 
situations. In this paper Ti, will be I 2 or M. The space of functions from H), 

whose generalized derivatives up to order m are also in is denoted by 

By definition = Lp^4R‘‘,n). The norm \f\w^^ of / in 

U) is defined by 

(1-1) I/IE.-, = E i^’”/r,„,. 

|a|<m 

where 11“ := for multi-indices a := {ai, ...,ad) € {O,!,...}"^ of length 

|a| := ai +a 2 ... -I- ad, and Dif is the generalized derivative of / with respect to 

a;® for i = 1, 2..., d. We also use the notation Dij = DiDj and Df = {Dif ,..., Ddf). 
When we talk about “derivatives up to order m” of a function for some nonnega¬ 
tive integer m, then we always include the zeroth-order derivative, i.e. the function 
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itself. Unless otherwise indicated at some places, the summation convention with 
respect to repeated integer valued indices is used throughout the paper. The con¬ 
stants in our estimates, usually denoted by N, may change from line to line in the 
calculations, but their dependencies will always be clear from the statements. 


2. Formulation and localization error estimate 


Consider the stochastic equation 

OO 

(2.1) dut{x) = {Ltut{x) + ft{x)) dt + ^{M^ut{x) + gt{x)) dw^ 

k^l 

on {t,x) € [0,T] X K.'^ =: Ht, with initial condition 

(2.2) uo{x) = ipix), X G 

Here / and g = are functions on U x Ht with values in K and I 2 , respectively, 

and L and are second order and first order differential operators of the form 

Lt = af{x)D^Dj+bl{x)D^ + Ct{x), = af [x)Di + g!l{x), fc=l,2,..., 

where the coefficients a*-l, 6*, c, and /i^ are real-valued functions on U x Ht for 
i,j = 1, 2, ...d, and integers A: > 1. 

For an integer m > 0, p G [2, 00 ), and d G R the following assumptions ensure 
the existence and uniqueness of a fU^^-valued solution (Mt(-))te[o,T] • 

Assumption 2.1. For P dt i^) dx-almost all (w, t, x) G Q x [0, T] x 

al\x)z’‘z^ > 0 

for all z G R'^, where 

= 2a*^' - 


This is a standard assumption in the theory of stochastic PDFs. Below we 
assume some smoothness on a in x, and on all the coefficients and free terms of 
problem (I2.ip - (j2.2() . We will also require that the nonnegative symmetric square 
root p := i/a possesses bounded second order derivatives in x. Concerning this 
assumption we remark that it is well-known from [3] that p is Lipschitz continuous 
in a; if q; is bounded and has bounded second order derivatives, but it is also known 
that the second order derivatives of p may not exist in the classical sense, even if a 
is smooth with bounded derivatives of arbitrary order. 

Assumption 2.2. (a) The derivatives in x G of up to order max(TO, 2) are 
V ® B{W^)-measurable functions, bounded by K for all i,j G {1, 2,..., d}. 

(b) The derivatives in x G R'^ ofP and c up to order m are -measurable 

functions, bounded by K for all i G {1,2, ...,d}. The functions ct* = and 

p, = ore l 2 -valued and their derivatives in x up to order m-\-l are 'P^BfMf')- 

measurable h-valued functions, bounded by K. 

(c) The derivatives in x G of p = y/a up to order m -G 1 are V 0 yB(R'^)- 
measurable functions, bounded by K. 


Assumption 2.3. The initial value, if is an Hq- measurable random variable with 
values in The free data, ft and gt = are predictable processes with 

values in and W^J'^{l 2 ), respectively, such that almost surely 

rT 


(2.3) 


KutAT) ■■= Aw; 


{\ft 


t\w" 


\9t 


w" 


n) dt < 


OO. 
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Definition 2.1. A 14^^ ^-valued function u, defined on [0, T] x fJ, is called a solution 
of (I2.1D - (I2.2D on [0,T] if u is predictable on [0,r] x Q, 

/ |ut|Li dt<oo{a.s.), 

Jo P’” 

and for each ip € for almost all a; € 

(ut.ip) ={'ip,p)+ f {-{alW^Us,Djip) + {blD^Us + CsUs +fs,(p)}ds 

Jo 

(2.4) + [ {a"DiUs +pIus+ gl,tp)dwl 

Jo 

for all t G [0,T], where b‘ = b‘ — Dja^^ , and (•, •) denotes the inner product in the 
Hilbert space of square integrable real-valued functions on 

The following theorem follows from Theorem 2.1 in [B]. 

Theorem 2.2. Let Assumvtions \2. 1[ lil.gl (a)-(b), and \2.S\ with m > 0 hold. Then 
there exists at most one solution on [0,T]. If Assumvtions lKH\2.2Y a}. and \2. A hold 
with m>l, then there exists a unique solution u = {ut)te[o,T] on [0,T]. Moreover, 
u is a -valued weakly continuous process, it is a strongly continuous process 

with values in and for every q > 0 and n € {0,1,..., m} 

(2.5) E sup \ut\l. <NEICl^^{T), 

te[o,T] 

where N is a constant depending only on K, T, d, m, p, d, and q. 

Remark 2.3. Theorem 2.1 in [B] covers the d = 0 case. We can reduce the case 
of general d to this by rewriting the equation for u as an equation for u(x) = 
u(x)(l-|- (see e.g. (3). It is easily seen that the coefficients of the resulting 

equation still satisfy the conditions of the theorem. 

Definition 2.4. A -measurable random field u on Ht is called a classical 

solution of if along with its derivatives in x up to order 2 it is continuous 

in {t,x) G Ht, it satisfies (I2.1I) - (I2.2I) almost surely for all {t,x) G Ht, and there 
exists a finite random variable ^ and a constant s such that almost surely 

\D°'ut{x)\ < .^(1 -b lx])'* for all (t,x) G Ht and for |a| < 2. 

Corollary 2.5. Let assumptions of Theorem 12.21 hold with m > 2 -b d/p. Then 
there exists a unique classical solution u to (CT)-(I^. 

Let us refer to the problem (I2.1I) - (I2.2I) as Eq(S)), where U stands for the “data” 

D = {ip,a,b,c,a,g.,f,g) 

with a = b = (R), a = (cr^*), g = {g^) and p. = (/r^). We are interested in 

the error when instead of Eq(S)) we solve Eq(S) with 

2) = {'4i,d,b,c,d,p,f,g). 

Assumption 2.4. Almost surely 

(2.6) D = £) on [0, T] x {x G : |x| < R}. 
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The main example to keep in mind is when each component of S) is a truncation 
of the corresponding component of S) (see Section [3] below). Let 

Br = {x € : \x\ < R} 

for R> 0. Define ^ ^{T) as JC^ ^ ^{T) with i/i, / and g in place of / and g, 
respectively. The main result reads as follows. 


Theorem 2.6. Let v € (0,1) and let Assumvtions \2.1\ lii.ill (h)-(c) and \2.d\ hold 
with m > 2 + d/p and "i? G R for S and S. Let also Assumption \2.4\ hold. Then 
Eqf'D) and Eq/Ti) have a unique elassieal solution u and u, respeetively, and for 
q > 0, r > 1 


(2.7) E sup sup \ut{x) — Ut{x)\'^ < Ne E^/’'(/C^ 

te[0,T]a:6B,,R 




(T) 




where N and 6 are positive constants, depending on K, d, T, q, r, tt, p, and v. 


First we collect some auxiliary results. The following lemma is a version of 
Kolmogorov’s continuity criterion, see Theorem 3.4 of [4]. 


Lemma 2.7. Let x{0) be a stochastic process parametrized by and continuous in 
0 G D G R^, where D is a direct product of lower dimensional closed balls. Then 
for all 0 < a < 1, q > 1, and s > pja, 


Esup \x{6)\‘^ < N{1 + |D|) 

9 


sup(E|x(6i)|«")i/® + sup 

9 9^9' 


/E|x(6l) -x(6»')|«® 
V \e - 


1/T 


where N = N{q,s,a,p), and \D\ is the volume of D. 


Lemma 2.8. Let (at)tg[o_T] (/?t)tG[o,T] be Et-adapted processes with values in 
R'^ and l 2 (R'^)) respectively, in magnitude bounded by a constant K. Then for the 
process 

(2.8) Xt = f Ois ds -\- 

Jo 

there exist constants e = e{K,T) > 0 and a N = N{K,T) such that 


1 

JO 


dUwl ie[o,r] 


Esupe^l^‘1'' < N. 

t<T 


Proof. A somewhat more general lemma is proved in [23] . For convenience of the 
reader we give the proof here. By Ito’s formula 

Yt := = 1 + / + 2a,A, 

Jo 


-\-2\f5sXg'\^ — /r|p} ds + mt 

for any /i £ R, where (wt)tg[o,T] is a local martingale starting from 0. By simple 
inequalities 


2aA + 2|/3A|2 < \a\^ + |Ap + 2\I3\^\X\^ < + {2K‘^ + l)|Ap. 

Hence for p = {2K^ + 1) and for a stopping time t < T we have 


EK 


tATn ^ 


< 1 + 2K^ 
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for Tn = T A Pm where is a localizing sequence of stopping times for m. 

Hence, by Gronwall’s lemma, 




< e 


2K-‘T 


where N is independent of n. Letting here n -A oo, by Fatou’s lemma we get 




2„-^T 




< e 


K^T 


for stopping times t < T. Hence applying Lemma 3.2 from [ 8 ] for r G 
obtain 


t<T 


— 1—r 


( 0 , 1 ) we 


□ 


To formulate our next lemma we consider the stochastic differential equation 

(2.9) dXs = asiXs) ds + P^iXs) dwt 

where a and /3 = (/?^) are V ® S(R'^)-measurable function on H x [0,T] x with 
values in R'^ and ^ 2 (R'^) such that they are bounded in magnitude by K and satisfy 
the Lipschitz condition in a; € R'^ with a Lipschitz constant M, uniformly in the 
other arguments. Then equation (12.91) with initial condition Xt = x has a unique 
solution X*’^ = {Xl’^)s^yt^T] for any t S [0,r] and x G R'^. 

Remark 2.9. It is well known from m that the solution of (12.81) can be chosen 
to be continuous in t,x,s. In the following, by Xl'^ we always understand such a 
continuous modification. 

Lemma 2.10. Set — x. There exists a constant S = 6{d, K, M, T) > 0 

such that for any R, 

(2.10) E sup sup < 7V(l + i?'^+i/2), 

0<t<s<T |a;|<R 

and for any R and r 

(2.11) P{ sup sup >r) <7Ve-'^’’'(l + i?‘^+i/2), 

0<t<s<T |x|<R 

where N = N{d,K, M,T). 

Proof. It is easy to see that (12.101) implies (12.lip , so we need only prove the former. 
For a fixed 6, to be chosen later, let us use the notations f{y) = and 7 = 

2{d + 2) + 1. By Lemma [2.71 we have 

E sup sup /(1‘’^) < 7V(1 + i?'^) sup sup (E/'>'(X*’“))i/'>' 

0<t<s<T |a;|<fi 0<t<s<T \x\<R 


(2.12) +7V(l + i?^) sup 


0<t<s<T |a:|<R \^{\t — t' \'^ + \s — s' \'^ + \x — x' \^)'1 / j 

0<t’<s'<T |x'|<R 

The first term above, by Lemma [2^ provided S < e/ 7 , can be estimated by NR'^. 
As for the second one, 

fixin - fixl^) = f dfidxr + (1 - dd. 

Jo 


E\f{xi’-) - f{ry)p 


1/7 
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Notice that |V/(?/)| < N{S)f^{y), therefore by Jensen’s inequality and Lemma [2^ 
again, provided S < we obtain 

Now the the right-hand side can be estimated by standard moment bounds for 
SDEs, see e.g. Corollary 2.5.5 in [15], from which we obtain 

-N ^.11 \ 1/ (27) 

(|t — t'P-I-|s — s'P-I-|a; — y “ 

□ 



Proof of Theorem \2.6\ 

Throughout the proof we will use the constant A = X{d,q), which stands for a 
power of i?, and, like N and d, may change from line to line. Clearly it suffices 
to prove Theorem 12.61 with e~^^ in place of e~^^ in the right-hand side of 
inequality (EJj). We also assume first that g > 1 and d = 0. 

The main idea of the proof is based on stochastic representation of solutions to 
linear stochastic PDEs of parabolic type, see [18], [19], and [21]. This representa¬ 
tion can be viewed as the generalization of the well-known Feynman-Kac formula 
and is derived as follows. First, we consider an equation which differs from the 
original only by an additional stochastic term driven by an independent Wiener 
process. The new equation is fully degenerate and taking conditional expectation 
with respect to the original filtration of its solution gives back u. On the other 
hand, the method of characteristics allows us to transform the fully degenerate 
equation to a much simpler one. This provides a formula for the representation of 
u, and, more importantly for our purposes, allows us to compare u and u on the 
level of characteristics. 

Recall that p = {pY{x))f^^i is the symmetric nonnegative square root of a = 
{2a^^ — and p is the symmetric nonnegative square root of d = {2dX^ — 

Then due to Assumption [Td] p = p almost surely for all t € [0, T] and 
for \x\ < R. Let {Wf)t>o,r=i...d be a d-dimensional Wiener process, also independent 
of the (T-algebra Pco generated by Tt for t > 0. Consider the problem 

dvt{x) ={LtVt{x) + ft{x)) dt + {M^vt{x) + gt[x)) dw'f 

(2.13) +J\ffvt{x)dwl 

(2.14) vo{x) ='ipix), 

where A/’*' = p'^^Di. Then by Corollary 12.51 (I2.13(l - (|2.14(1 has a unique classical 
solution u, and for each t £ [0, T] and x almost surely 

(2.15) Mt(a;) = E(ut(a;)|J't). 

Together with (12.131) let us consider the stochastic differential equation 

(2.16) dYt = l3t{Yt)dt-a^{Yt)dw^-pl{Yt)dwl 0 < t < T, Yo = y, 
where 

Ptiy) = -bt{y) + af{y)Dia^{y)+pl\y)D,pl{y) + a^{y)p.^{y), t £ [0,T], y £ 
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and (T^, p’' stand for the column vectors (cr^^,..., , p'^’'), respectively. 

By the Ito-Wentzell formula from [16], for 

Ut{y) := Vt{Yt{y)) 

we have (to ease the notation we omit the parameter y in Yt{y)) 

dvt{Yt) = {LtVtiYt) + ft{Yt)) dt + {M^Vt{Yt) + g^{Yt)) dw^ +Af[vtiYt) dw^ 
+[l3lD,VtiYt) + al^D,^Vt{Yt)) dt - al’^DMYt) dw^ - M[v{Yt) dw^ 

(2.17) - <jfD,{M^vt{Yt) + g>^{Yt)) dt - J\f[J\r[vt{Yt) dt. 

Due to cancellations on the right-hand side of (I2.17P we obtain 

dUt{y) ={^t{Yt{y))Ut{y) + MYtiv))} dt 

+ {t^t(Yt{y))Ut{y) + gt{Yt{y))}dw’^, C/o(y) = tp{y), 

where 

7 t(x) := ct{x) - erf (a;)APt^(a;), = Mx) - a^\x)D,g^. 

Notice that in the special case when / = 0, p = 0, c = 0, p = 0 and '(/’(a;) = a:* for 
i G {1, ..., d}, we get vl{Yt{y)) = y® for i = 1, ..., d, where u® is the solution of (12.1311 - 
(I2.14P with f = c = 0,g = g = 0,a = 0 and ’tp{x) = x^. Hence for each t G [0, T] the 
mapping y —>■ Yt(y) G has an inverse, Yf~^, for almost every w, and the mapping 
X —>■ Vt(x) = (Vf(x))f^i, defined by the continuous random field (vt)(t,x)eHT gives 
a continuous modification of Yf~^. Also, we can write Vt(x) = Ut(Yf.~^). 

As we shall see, due to the data being the same on a large ball, the character¬ 
istics Y and Y agree on an event of large probability. This fact and the above 
representation will yield the estimate (|2.7|) . Set Ut(y) = vt(Yt(y)), where vt(x) and 
Yt(y) are defined as vt(x) and Yt(y) in p.l3p - ()2.14p and (12.1611 . respectively, with 
S) and p in place of S) and p. 

Introduce the notations B/j = [0, T] x Bji and = B/j fl Since u and u 

are continuous in both variables, 

(2.18) sup \ut{x) - ut{x)\ = sup \ut{x) - ut{x)\ 

(t,x)6B„H (t,x)GA„_R 

Let !/' = (! + i')/2 and define the event 


H := 

sup \Y^ ^(x) > v'R 

U 

sup |yt(x)|>i? 


(t,x)gB„H 


(t,a:)eB„,R 


Then 

= [rrfx) G B,,R,\/it,x) G B,r] n [Ytix) G BR,\fit,x) G B,>r] , 
and thus on H‘^ 

Yt{x) = Ytix) for {t,x) G B^>r, 

Yf'^ix) = Y~^{x) for {t,x) G B„r , 

and consequently, 

vtix) = vt(x) for (t,x) G B,,r. 

Therefore, by (12.1511 and (I2.18L and by Doob’s, Holder’s, and the conditional Jensen 
inequalities, 

E sup \ut{x) — utix)]^^ < E sup |E(l/i- sup \vTix) — VTix)\\Bt)\‘^ 
(t,a:)GB„H te[0,T]nQ {t,x)^A^r 
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(2.19) 

< \iPiH))^R'E^Ri sup \vrix) vAx)r') 

(2.20) 

< ^Ap{h)Y/^Vt 

q-l 

with 

sup \Vr{x)A YE^/^ sup \Vrix)A 
{t,x)^Ht {t,x)^Ht 

for r > 1, r' 

= r/(r — 1 ), provided g > 1. By Theoreml2.2l 

(2.21) 

Et < NE^/^iKAj,AT)+^'‘r::,pAT))- 


We can estimate P{H) as follows. Clearly, 

P{H)<P{ sup \Y^^{x)\> v'R) +P{ sup \Yt{x)\> R) =: Ji +J 2 - 

For Yt{x) = Yt(x) — cc by (12.111) we have 

J 2 < Pi sup |yt(a:)| > (1 - i^')R) < 

Also, we have 

00 

Ji < ^ P(3(t, x) e [0, T] X iB 2 i+i^>R \ B 21 v>r) : \Ytix)\ < vR) 

1=0 

CO 

<^P( sup \Ytix)\>i2‘R'-R)R). 

Using f|2.1ip again gives 

00 

1=0 

We can conclude that 

( 2 . 22 ) P{H)<Ne-^^\ 

where N and 6 are positive constants, depending only on d, K and T. 

Combining this with (12.201) and (12.211) we can finish the proof of the theorem 
under the additional conditions. 

For general d one applies the same arguments as in Remark 1 2.3 1 Finally (12.71) for 
the case q € (0,1] follows easily from standard arguments using Lemma 3.2 from 

0 . □ 


3. An APPLICATION - FINITE DIFFERENCES 

In this section we apply Theorem 12.61 to present a numerical scheme approxi¬ 
mating the initial value problem <STT\i-IST^. We make use of the results of [7] on 
the rate and acceleration of finite difference approximations in the spatial variable, 
which, together with a time discretization and a truncation - whose error can be 
estimated using Theorem l2.6l - yields a fully implementable scheme. We shall carry 
out the steps of approximation in the following order: spatial discretization by fi¬ 
nite differences, localization of the finite difference scheme, and discretization in 
time via implicit Euler’s method. This of course requires an analysis of the Euler 
scheme, to present an error estimate for it, which does not depend on the spatial 
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mesh size and the localization. Furthermore, in our full discretization scheme we 
shall incorporate Richardson’s extrapolation, which will allow us to improve the 
accuracy of the scheme in the spatial mesh size h. 

First we introduce the finite difference approximation in the spatial variable for 
(lOl-dO). To this end, let Ai C be a finite set, containing the zero vector, 
satisfying the following natural condition: Aq := Ai \ { 0 } is not empty, and if a 
subset A' C Aq is linearly dependent, then it is linearly dependent over the rationals. 
Let h > 0, and define the grid 

n 

= {h 'y ] Ai : Ai S Ai U — Ai, n = 1,2,...}. 

i=l 

Due to the assumption on Ai, Gh has only finitely many points in every ball around 
the origin in Define for A € Aq U — Aq, the finite difference operators 

+ hX) - 2(f(x) + <f(x - hX)), 

and let Sq stand for the identity operator. To approximate the Cauchy problem 
( 1 ^ - d^ . for h > 0 we consider the equation 

OO 

(3.1) dut{x) = (L’^utix) + ft{x)) dt + '^{M^’^utix) + g^{x))dw'^ 

k=l 

on [0,T] X Gfi, with initial condition 

(3.2) uo(x) = ip{x), 

where and difference operators of the form 

L’lix)= Mf’'=(x) = ^ b,^''=(:r)d^ fc = l,2,..., 

A,k;^A.i AgA.1 

with some real-valued V 0 S(]R^)-measurable and on D x [0,T], such that 

(3.3) |a^''(a:)| < A" and ^ |b^'=(a;)p < AT^ 

k 

for all A,K G Ai, t G [0,T], a; G K'’* and uj G XI, where AT is a constant. 

Remark 3.1. Flere "0, / and g are the same as in (I2I1)-(E21) and we will assume that 
they satisfy Assumption 12.31 with m > d/2, p = 2 and d = 0. Thus by Sobolev’s 
embedding of into Cb, the space of bounded continuous functions, for all uj we 
can find a continuous function of x which is equal to ip almost everywhere, and for 
each t and u! we have continuous functions of x which coincide with ft and gt for 
almost every a; G Here and in the following we always take such continuous 
modifications if they exist, thus we always assume that ip, ft, and gt are continuous 
in X for all t (for g = {g^)'^i this means, as usual, continuity as a function with 
values in 12 ). In particular, terms like ft{x) in (13.11) make sense. We note that for 
TO > d/2 one can use Sobolev’s theorem on embedding to Cb to show also that 
if Assumption 12.31 holds with m > d/2, p = 2 6 = 0, then 

^ \ip{x)\^h<^+ r \Mx)\^h'^+ E 

xGGh ^ x^Gfi x^Gh fc 


LOCALIZATION ERRORS 


11 


< NU{x)\\l^ +N I ||/t||^ + ^||gf||^dt <oo(a.s.), 

•^0 k 

with a constant N = N{Ao,d), where || • ||m := | • |m, 2 .o- (See Lemma 4.2 in [T0].i 

Clearly, for ip G and A 0 

5^ip[x) —>■ d\ip{x) := X^Diip{x) as h —>• 0. 

Thus, in order to approximate L and by and respectively, we need 

the following compatibility condition. 

Assumption 3.1. For every i,j = l,...,d, fc = 1,... and for P dt (g) dx-almost 
all {uj,t,x) G n X [0,T] X 

a-^ = > a A , b = y. )a , c = a , 

A,k.gAq AgAq 

cr*'= = / = b°’'=. 

A^Aq 

For each x G G/i equation (EH) is a stochastic differential equation (SDE), i.e., 
in general, (I3.1I) - (I3.2I) is an infinite system of SDEs. To replace this with a finite 
system we make the coefficients, together with the free and initial data, vanish 
outside of a large ball by multiplying them with a cutoff function C/j, which satishes 
the following condition. 

Assumption 3.2. For an integer m > 0 and a real number R > 0 the function 
is a continuous function with compact support on W^, such that C,{x) = 1 for |a;| < R 
and the derivatives of up to order m+1 are continuous functions, bounded by a 
constant C. 

In this way we replace (I3.1I) - (I3.2I) with the system of SDEs 


(3.4) dut{x) = {LYutA) +ft^i 

x))dt+ {Ml/’^'’'ut{x) + g^’'"{x)) dw /, 

tG [o,r]. 

with initial condition 





(3.5) 



uo{x) 

= i}^{x), 


for X G 

Gh DsuppCr, 

where supp Cr is 

the support of Cr, 


(3.6) 



R,k\ 

9t )■ 



and 







lY-.= E 


, 

,R,k _ ^x,R,k^h^ fc = l,2. 

. . . , 





AeAi 


with 






(3.7) 


^Xk,R . 

= C|a' 

for A, K G Aq, 


(3.8) 

(a0/..fl, ^xo,R 

,f,A.R.fc) _ 

^ (CRa° 

'',(Ra^°,CRb^^^) forA,KGAi, 

fc > 1. 


At this point our approximation is a hnite dimensional (affine) linear SDE, whose 
coefficients are bounded by K owing to (j3.3p . and furthermore, by virtue of Remark 
13.11 for each x, f^{x) and g^{x) are square integrable in time under Assumption 
13.21 and 13.51 below with m > d/2. 
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Hence ()3.4ll - (j3.5l) has a unique solution 

{ui'^{x) : X € G/, nsuppCR}tG[o,T], 

by virtue of a well-known theorem of Ito on finite dimensional SDEs with Lipschitz 
continuous coefficients. The approximation of such equations are well studied, 
various time-discretization methods can be used, each of them with their own ad¬ 
vantages and disadvantages. Here we chose the implicit Euler method, formulated 
as follows. 

We take a mesh-size t = T/n for an integer n > 1, and approximate (I3.1I) - (I3.2I) 
by the equations 

(3.9) wo(a;) = ip^{x), 

Vi{x) = Vi-i{x) + {L^fviix) + 

OO 

(3.10) i = l,2,...,n, 

k^l 

for a: € G/, U suppi^R, where 

Remark 3.2. In many applications, including the Zakai equation for nonlinear fil¬ 
tering, the driving noise is finite dimensional. If this is not the case, one needs 
another level of approximation, at which the infinite sum in (13.101) is replaced by 
its first m terms. We shall not discuss this here. 

Remark 3.3. As mentioned before, Euler approximations for SDEs are very well 
studied. Therefore, while it is far from immediate that the error is of the desired 
order, independently of h and R, the implementation of the scheme goes as usual, 
see e.g. [14] and its references. 

To prove to solvability of the fully discretized equation ()3.9|) - (|3.10ll and estimate 
its error from the true solution of (I2.1I) - (I2.2D on the space-time grid we pose the 
following assumptions. As for the following we confine ourselves to the L 2 -scale, 
without weights, we use the shorthand notation || • ||„ = | • \ w^g, II ' II = II ’ llo- 

Assumption 3.3. For all {uj^t,x) £ D x [0,T] x 

( 2 ^^” - > 0 

A,k^A.o 

for all z = (z'^)AeAo, G R. 

In the following assumptions m and I are nonnegative integers, as before, and 
will be more specified in the theorems below. 

Assumption 3.4. The derivatives in x ofa^'^ up to order max(rn, 2) are P0H(M‘^)- 
measurable functions, bounded by K for all X,k G Ai. The derivatives in x of 
b^ = (b^’’))lT]^ up to order m -\- 1 are V ® -measurable l 2 -valued funetions, 

bounded by K, for all A G Ai. 

Assumption 3.5. The initial value, if is an ipQ-measurable random variable with 
values in W™. The free data, ft and gt = (ff^)^i are predictable processes with 
values in IE™ and respectively, such that almost surely 
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Assumption 3.6. There exists a constant H such that 

E||/t - M\f+E\\gs - gtWl, < H\t - s\, nftWhi + E||<?*ll ?+2 < H 
for all s^t G [0, T], and 

- af (x))p < H\t - s\, ’'(^) " < H\t „ 

k 

for all w G n, X G , s,t G [0,T] and multi-indices a and /3 with |a| < I and 

\f3\<l + l. 

Remark 3.4. If Assumptions I3.1| and 13.31 hold then 

(2a*^' - k^)z^ z^ 

A,k^A.q AgA.0 /^GA.q 


A, K 

for all z = {z^,...,z‘^) G i.e., Assumption 12.11 also holds. Clearly, Assumptions 
13.II and 13.41 imply Assumption 12.21 (al-lbl. Notice that Assumption 13.51 is the same 
as Assumption 12.31 with p = 2 and i? = 0. Thus if Assumptions 12.21 (c), and 13.11 
through [331 hold with m > 2 -h d/2, then by virtue of Corollary 12.51 equation p.lll 
with initial condition (12.21) has a unique classical solution 

u = {ut{x) : t G [0,r],x G M'^}. 

Now we are in the position to formulate the first main theorem of this section, 
with the notation G// = Gh H Bji. 

Theorem 3.5. Let I > d/2 be an integer. Let Assumvtions lKT] throuqh \3.4\ hold 
withm > 4 + /, and let Assumvtions \2.S\ (c) and Assumvtion \3.6\ hold withm > 2-\-l 
and with / + 1, respectively. Then ifr is sufficiently small, then for any h > 0, R > 1 
the system of equations (13.91) - (13.101) has a unique solution Moreover, 

for any v G (0,1), g > 1 we have 

E max max lurAx) — 
i=0,....na:eG^« 


(3.12) < -h N2ih^ + t)(1 E/C^), 

with constants Ni and 5 > 0 depending only on K, d, T, C q, v and Ag, and a 
constant N 2 = N 2 {K, T, d, C, H,Aq). 

As mentioned above we want to have approximations with higher order accuracy 
in h by extrapolating from recall the method of Richardson’s extrap¬ 

olation. This technique, first introduced in [22], allows one to accelerate the rate 
of convergence by appropriately mixing approximations with different mesh sizes, 
given that a power expansion of the error in terms of the mesh sizes is available. 
We shall use this, based on results of |7], to obtain higher order approximations 
with respect to the spatial mesh size h. To formulate the extrapolation, let r > 0, 
V be the (r -I- 1) x (r -|- 1) Vandermonde matrix Rb = 

(3.13) (co,ci,...,c,):=(l,0,...,0)V-\ 
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and define 

(3.14) 

i^O 

where denotes the solution of ()3.9I1 - (I3.10I1 with h/2® in place of h. As we 

shall see, even by mixing only two approximations with different mesh sizes, that 
is, setting r = 1, the extrapolation increases the order of accuracy in h from 2 to 4. 
The second main result of this section is the following. 

Theorem 3.6. In additions to the assumptions of Theorem \3.5\ let Assumptions 

I,*?..21 \3.4\ and \3.5\ hold with m > 4r + 4 + Z. Then for the extrapolation 

have 

E max max \uri{x) — {x)\^ 

(3.15) < + iV2(/i2(2r+2) 

for any v £ (0,1) and q > 1, with constants Ni and (5 > 0, depending only on K, 
d, T, C, V, q and Aq, and a constant N 2 = N 2 {K,T,d,C, H,r, Aq). 

These theorems will be proved by using Theorem l2.6l some results from [7], sum¬ 
marized below in Theorem 13.91 and the error estimate for the time-discretization, 
established in Theorem 13.101 below. 

Example 3.7. Consider the equation 

dut{x) = siv?{x)D‘^ut{x)dt -|- sm{x)Dut{x)dwt 

for {t,x) € [0,1] X R, where (wOtefo,!] is a 1-dimensional Wiener process, with the 
initial condition 

Uo{x) = (1 -I- x'^y^. 

The choice of localizing function C/j is quite arbitrary, for the sake of concreteness 
we take Cr{x) := f(x -I- 2 -|- i?) — f{x — 2 — i?), where 

2 2 

f{x) := — arctane'^/i^”® i for |a;| < 1, and f{x) := l[i^oo)(2;) for |a;| > 1, 

while noting that in practice a simple mollified indicator of [—R,R] may be more 
favourable. Notice that Cr{x) = 1 for \x\ < R and supp = [—3 — R, 3 -I- R]. 

For an integers j > 1 and n > 1 we set h = i?/(10j) and r = 1/n. To use the 
extrapolation with r = 1 in (13.141) . we need to solve two discrete equations with 
spatial mesh sizes h = h,/i/2, and mix them according to (13.141) . where one can 
check that the coefficients are cq = —l/3,ci = 4/3. Following the steps outlined 
above, the discrete equation we arrive at is 

uf'^''^{kh) = a{kh){S^S^uf’^){kh)T + b{kh){S^ufz^){kh){wri — Wr(i-i)) 

for i = 1, 2,..., n and fe = 0, ±1,..., ± [(3 -I- i?)/h], with the initial values 

UQ’^’'^{kh) = (1 + k^yy^ (R{kh), 

where a(a;) = Cfl(®) sin^(a;), b(a;) = Cr(®) sin(a;), and 

S^(j){x) = {2hy^[(j>{x + h) — 2(j){x) + (j){x — h)]. 

For each h one can solve the above equation recursively in i. 
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Taking the resulting solutions and and setting 

^R,h,r ^ (4/3)y«.V2.r _ 

we can conclude that the error Emaxi \uiT{kh) — {kh)\^^ where i runs over 

0,1,..., n and k runs over Oil,..., ±0.9R/h, is of order i h® i r. 

Before summarising some results from [7] on finite difference operators L^, 
and stochastic finite difference schemes in spatial variables we need to make an 
important remark. 

Remark 3.8. The concept of a solution of (I3.ip - ()3.2I1 . as a process with values in 
h.h^ the space of functions (/) : G/j —?► K with finite norm ^ = J2xeGh 
is straightforward. One can, however, also consider on the whole space, 

that is, for {t,x) € Ht- Namely, when Assumptions 13.41 and 13.51 hold, then we can 
look for an J't-adapted L 2 -valued solution (u^)te[o,T] such that almost surely for 
every t € [0, T] 

pt pt 

(3.16) u1=i;+ + f,)ds + J2 M 

Jo 

in the Hilbert space L 2 , where the first integral is understood as Bochner integral 
of L 2 -valued functions, the stochastic integrals are understood as Ito integrals of 
L 2 -valued processes, and the convergence of their infinite sum is understood in 
probability, uniformly in t G Thus, by a well-known theorem on SDEs in 

Hilbert spaces with Lipschitz continuous coefficients, equation (13.161) has a unique 
L 2 -valued continuous J't-adapted solution = (uJ‘)tg[o,T]- We refer to such a 
solution as an L 2 -valued solution to (IXTD-(IO). We can view equation (I3.16|) also 
as an SDE in the Hilbert space IT™, and by the same theorem on existence and 
uniqueness of solution to SDEs in Hilbert spaces, we get a unique lT™-valued 
continuous Jr-adapted solution to it. Consequently, if Assumptions 13.41 and 13.51 
hold with m > 0 , then u = (ut‘)te[o,T]i the L 2 -valued solution to (iTTD-dT^ is a 
lT™-valued continuous J't-adapted process. Also, it is straightforward to see that 
these two concepts of solutions are “compatible” in the sense that if m > d/2, then 
the restriction of the L 2 -valued solution to Gh solves (EH)-( 021 ) as an / 2 ,ii-valued 
process. Note that (I3.4I) - (I3.5I) is a special case of the class of equations of the form 
EH-EH, and so the above discussion applies to it as well. 

The analogous concepts will be used for solutions of (13.911 - (13.101) . Namely, a 
sequence of L 2 -valued random variables is called an L 2 -valued solution 

to (j3.9l) - (l3.10ll if ufis J^ir-measurable for i = 0 , 1 , 2 , ..., n, and the equalities hold 
for almost all x € for almost all ui € fl. If (u/’ )(hQ is an L 2 -valued solution to 
(I3.9I) - (I3.1QI1 such that € IT™ (a.s.) for m > d/2 for each i, then it is easy to 
see that the restriction of the continuous version of to Gh H supp Cr for each i 
gives a solution {vi(x) : x € Gh nsuppCR,i = 0,1,...,n} to (I3.9D - (I3.10D . 

Theorem 3.9. Let Assumvtions ld.31 \3.4\ and \3.5\ hold with an integer m > 0. Then 
(a) For any </> G IT™ and I 7 I < m 

2(D'^</>, ^ < NUWl 

k 

for all oj ^ Vt and t G [0, T] with a constant N = N{K^ m, Aq); 
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(b) There is a unique L 2 -valued solution of (I3.1ll - (j3.2|) . It is a Wip-valued 
process with probability one, and for any q > 0 

Esup||iif||^<iVE/C^ 

t<T 

with a constant N = N{K, m, d, Aq, q, T); 

(c) If for an integer r > 0 Assumvtions fXTl through \ 3.5\ with m > 4r + 4 + d/2 
hold, then for any q > 0 

Esup max \ut{x) - 

tKT^eGh 

where u is the classical solution to (I2.1I) - (I2.2I) . N is a 

constant depending only on K, d, T, q, Aq and m. 

As discussed above, under Assumptions 13.41 and I, 3. 51 we have a unique L2-valued 
solution = [u^’ ^)te[o.T] to dSa-dSA]), and it is a continuous bF™-valued pro¬ 
cess. 


Theorem 3.10. (i) Let Assumvtions \S.‘A through \3.5\ hold with m > 0. Then for 
sufficiently small r there exists for all h and R > 0 a unique L 2 -valued solution 
yh,R,T (|3.9p - (|3.l0p such that € IT™ for every z = 0,1, ...,n and uj £ LI. 

(ii) If Assumption\KM holds with some integer I > 0 and Assumptions\3fl\ through 
\3.5\ hold with m = I Z, then 

(3.17) maxEllu^f - < Nt{1 + E/C^), 

i<n 

with a constant N = N{K, C, H, d,T, I, Aq). 

(Hi) Let I > 0 be an integer. If Assumvtionholds with I 1 in place of I, and 
Assumvtions \3J\ throuah \3.5\ hold with m = I A, then 

(3.18) EmaxWuff-v^’^’^Wf < iVr(l+E/C^) 

i<n 

with a constant N = N[K, C, H, d,T, I, Aq). 


Proof. To prove solvability of the system of equations (l3.9() - (j3.10() we rewrite (13.1011 
in the form 


(3.19) 

{I - TL!lf)Vi = Vi^i 


OO 


k=l 


i = 1,2,... ,n, 


where I denotes the identity operator. We are going to show by induction on j < n 
that for sufficiently small r for each j there is a sequence of lT™-valued random 
variables {vi)l^Q, such that u* is JiT--measurable, vq = tp^ and p.l9l) holds for 
1 ^ ^ J- For j = 0 there is nothing to prove. Let j > 1 and assume that our 

statement holds for j — 1. Consider the equation 


(3.20) Du = X, 

where 

OO 

D := / - rLff, X := v,., + rU^-i) + 
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In the following we take k to be either 0 or m. It easy to see that D is a bounded 
linear operator from W 2 into W^, for each w G SI and r. Let us define the norm 
I • u in VL2" by 


d 

7 :| 7 |</« i—1 

where is a positive integer for each multi-index 7, jyj < k. Thus | • is a Hilbert 
norm which is equivalent to to || • ||k. We denote the corresponding inner product 
in by (, )«■ By virtue of Theorem 13.91 fal for all a; € SI and r we have 

(v, Dv)^ = |uU - v)^ > \v\l - tN\v\1, for all v G 

where the dependence of N is as in the theorem, in particular, it is independent of 
h,R. Consequently, for r < 1/A^ we have 

(u,]D>u)k > (5|u|^ for all v G W 2 , w G SI, 

where (5=1 — tN > 0. Hence by the Lax-Milgram lemma for every w G SI there is 
a unique v = Vk. £ W 2 such that 

(Dw^, (/?)« = {X, ip)^ for all (p G 

Since {Y, ip)^ = (Y, (/ — AYp>) for all Y G W 2 and ip G we have 

(Du^, (/ - A)''(/>) = {X, (/ - A)'‘(p) for all ip G C'o°°(K'^). 

Hence, taking into account that {(/ — A)™(^) : ip G C'“(R'^)} is dense in IT® = ^2, 
we get that Vm solves p.20l) in L 2 as well, so by uniqueness, Vm = vq. This means 
(13.201) has a unique solution v £ L 2 for every a; G SI, and v £ IT™ for every 
a; G SI. Since X and ]D>(/) are IT™-valued J^v-measurable random variables for every 
ip £ IT™, the unique solution v G IT™ to p.20l) is also J^T-measurable. This 
finishes the induction, and the proof of statement (i) of the theorem. 

For parts (ii) and (iii) < 00 may and will be assumed. Hence by Theorem 

[3H(b), 

(3.21) E sup < iVE/Cg < 00 for p = 0,1, 2,..., m. 

te[o,T] 

As Assumption 13.61 also holds, we have 

(3.22) E||V>||f+^ -HmaxE||/t,||^+^ -Hinax^E|cit';||^+,, < 00 

k 

with K = 0 in part (ii) and /c = 1 in part (iii), and hence, with the same k, 

(3.23) E||z;'‘’-^’^||f+^ < 00 for every I = 0,1,..., n. 

To start the proof of (ii), let us fix a multi-index 7 with length I7I =: g < 1. From 
equations (I3.9D - (I3.10D and (I3.1I) - (I3.2I) . we get that the error et := is 
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a 14^™-valued J>i-nieasurable random variable, i = 0 ,..., n, and is the L 2 - 

valued solution of the equation 


=Z?T'e,_i 


+ D^L^^ferT + 


D^Fs ds 


/r(i-l) 


+ + r D^G’: dw, 

J r(i—l) 


T(i-l) 

for i = 1 ,..., n, with zero initial condition, where 

F, := + /« - «2(t) := := (NJ + l)/n, 


G k -\jh.k h.R Ti/rn.ti.K n.FC , tc.K fc /i\ n/i\ \ ii / 

, ■=Mf UF +5t Ki(t) := (t) := [ntj/n. 


To ease notation we set 


fTl pTl 

&r-=L’lf, m'y-.= , d, -.= D'^Fsds, ©. :=/ D^G 

J J rd — l) 


dw^ 


/r(i-l) 


and by using the simple identity ||6|p — ||a|p = 2(6, 6— a) — ||6— a|p with 6 := D~^ei 
and a := D'^Ci-i, we get 

= 2{D^e,, + S,) + 2(i^'>'e„ + 6,) - 

= 2{D'^e^, D^^2,e^T + S,) + 2{D^e,-i,D^Tlt^e^-i^^ + 0.) 

+ 2(7^'^e, - + &^) - \\D^e, - 

= 2{D'^e^, D^&,ie^T + + 2{D^e^-l, + &^) 

+ + ©.f - 

(3.24) 

r(l) _|_ r(2) _|_ r(3) , r(4) ^(5) . r(6) 

— -^7 ~i” -t,• ~r J,• “r -i,• I -t \ -La 


with 


:=2T{D^ei,D^il,ei) 

/f) :=2{D'^e,,d,) 

/f) :=2(i4^e,_i,Z4^1Htie,_ief + ©*) 

/f) 

j(5) :=2(i4^1Htie._ief,©0 

4®^ :=ll®.f ■ 


By the Young and Jensen inequalities, and basic properties of stochastic Ito inte¬ 
grals we have 

(3.25) <T||7J'^e,f+ r-i||J,f <r||ZJ^e,f + / \\Fs\\lds, E/f^ = 0, 

J t( 2—1) 
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(3.26) 




l||^ 


pTi _ pTl _ 

(3.27) E/f)=E/ ^IIG^-fds. 

Jrir-l) V V 


^(*-1) ■ k 

By Ito’s identity for stochastic integrals 

E/f^=2E/ [ D'^m';:_iei-i{x)D'^G^,{x)dxds. 

J r (z—1) 

Here by integration by parts we drop one derivative from on the term 

D"^Gg, and then by the Cauchy-Schwarz-Bunyakovsky and Young inequalities we 
get 

Ejf) < T7VE||e,_i||2 + NE T ^ \\GX+i ds. 

Using (13.261) . by Theorem 13.91 (al we have 

(3.28) E/W+E/f) < iVTE||e,||2. 

Therefore, by taking expectations and summing up (13.241) over i from 1 to j < n, 
and over 7 for jyj < I, we get 

(3.29) E\\eX<Norj2n\er\\f + NoE i\\FX +J2\\^’'s\\hi) ds 

i=i do fc 

for j = where Nq = No{K,G,Ao,l,d) is a constant. Notice that due to 

(IT^ and (Irni) 

E||ei||f < 00 i = l,2, ...,n, 

and due to (|3.21l) and Assumptions 13. 4| and 13.61 we have 

i-T 


E 


f ||F,||2+^||Gj||2+,ds<oo. 
do 


Hence the right-hand side of inequality (13.291) is finite. Thus when r < 1/iVo, from 
(I3.29|) by discrete Gronwall’s lemma it follows that 

(3.30) E||e,||? < iVo(l - Noxym {\\FX + ^ ||G^|?+i) dt 

do k 

for j = 1,..., n. Now we are going to show that 

(3.31) EllA^llf+ ^E||G,^||?+i <iVr(E/Cf+3 + l) 

k 

for all t G [0,T] with a constant N = N{K, G,c,d,T, I, Ao). To estimate E||Ti||^, 
first notice that due to Assumption 13.61 

(3.32) E\\X-fX)\\l<NT, 
and due to Assumptions 13.2113.41 and 13.61 

-h,R h,R ph,R h,R \\2 


\LXuX - KfitKXf < ^Mt) + 2A,it) 


with 


Mt) := UlX - LXXtl)\\i ^ 


h,R 112 
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A2[t) := 


h,R/ h,R 


— U 


K2(t)' 


l?<^lh 


'K.2{t) 


\\l+2- 


By virtue of Theorem (b) and Assumption 13.2 

(3.33) 


EAi{t)<NT sup E\\ut’^\\f _^_2 < NTEK.f+ 2 - 

te[o,T] 


To estimate A 2 we show that 


(3.34) 


h,R 


-«^''ll?+ 2 <^l^-s|(l + E/C?+ 3 ) foralls,te [0,T]. 


To this end we mollify in the spatial variable. We take a nonnegative </> S 
C'“(R'^) supported on the unit ball in such that it has unit integral, and for 
functions (/?, which are locally integrable in a; € R'^, we define by 

■.= e~'^ [ (p{y)cj}{{x — y)/e)dy, x G R"^, for e > 0. 

Jr 

We will make use of the following known and easily verifiable properties of 
For integers r > 0 and e > 0, 

(3.35) - Lp\\r < £:||i^||r+i for (p G W 2 ~^\ 
and 

(3.36) <11^11., ||(A¥^)(^)||. = ||A<^(^)||. <^11^11., v^GWJ 

for i = 1,2, where W depends only on the sup and Li norms of -Dicf). Thus 
by (13.351) 

= Wut'"' - + E 

< +E ii?+i 

+EiiA«E- (Au*'‘’^)(^)|i?+i + Ell(^*“^''- 


2=1 


2=1 


||;+i 


h,Ru2 




h,R 
U. — U 


h,R 




2=1 


(3.37) +Ne\\\u’:-%+2 + \\u>:’%+2). 

Since satisfies (l3T])-(IO). for 0 < s < t < T we have 


(3.38) 

with 


h,R 


- w^^ll?+l < 2Bi + 2B: 


Bi := 


+ f^) dr 


i+i 


B2 := 


and using (I3.36P we have 


[ + g^'^) dwl 

J S 


/ + 1 


N. 


E \m4’^ - \\i, < ^B, + S 3 , 


2 = 1 


(3.39) 
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with 


^3 := X! II + A5.^’'=) dwl 

i=i 


i+i 


It is easy to see that 

Ei?! <{t-s) [ E||L^X"’" + /, 

J S 

(3.40) 

Using Ito’s identity, for B 2 we obtain 

Ei ?2 =E f J2 + 9r\\hi dr 

J S 7. 




< IV(t- s)^(supE||ut ||?+3 +supE||/t||f+i). 

t<T t<T 


(3.41) 

and in the same way, for B 3 we have 


< N{t - s)(supE||M^'^||^+2 + supE||5t||f+i), 

t<T t<T 


(3.42) 


EB 3 < N{t - s)(supE||u^'”||^+3 + supE||cit||f+2). 

t<T t<T 


From (|3.37p through (|3.4ip we get 

- Us’^f +2 < + {i- s)^e"^ + \t- s\)J 

for every e > 0, where 

J := supiEllu^-^ll^+g +E||5t||f+2 + E||/t||f+J. 


t<T 


Choosing here e := |t — gives 

Hence we obtain (|3.34p . since by Theorem 13.91 and Assumption 13.61 we have 

J < iVE/Cf+3 + H. 

From (13.341) we get the estimate 

EAi < 7Vr(E/C?3_3 + 1), 

which together with (13.331) and (13.321) shows that E||Ft||^ is estimated by the right- 
hand side of the inequality (13.311) . Similarly, by making use of (I3.34p . we can get 
the same estimate for E||Gt||^_|_ 3 , which finishes the proof of (13.311) . From (13.301) 
and (13.311) we have 

maxE||e,||2 < tN{1 - 7Vor)-"(E/Ce+3 + 1) 

j<n 

for T < Nq^, with a constant N = N{K, G, H,d,l, T, Aq). Hence noticing that 

lim (1 - Nqt)-^ = e^°^, 

n—>-oo 

we obtain (13.171) . which completes the proof of (ii). 

To prove (hi) notice hrst that by using (ii) with / -)- 1 in place of I we have 


E^||e,||2+i<A(l + E/C^) 


(3.43) 
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with m = I + A. Fixing a multi-index 7 as in (ii), with I 7 I = £< < I, we revisit (13.2411 
and observe that and by Theorem 13.91 

4 '^ < NrWe^Wl 

Thus summing up (|3.24l) over i = 1,..., j, and recalling also (I3.25P and (I3.26p . we 
obtain 

n p'j' 

E max WD'^aW^ < NtES" WaWf + N / EHi^.Hfds 
i<3<n Jo 


+N rE||G«||2ds + rf^E^||i4^antie,_if+E max^^/f). 

•^0 i=l fc -■ 

Hence, noticing that 

00 

J2\\D^^lle^-lr < N\\eJ\l„ 

k^l 

by using (|3.43p and (|3.3ip . we get 

3 

(3.44) E max ||i^^ei||^ < A^r(E/C^ + 1) +E max 

l<i<n 

2 = 1 

Clearly, 

^24"^ = r Z!jdwl j = l,2,...,n, 

i=i -^0 

where (^t')tG[o,T] is defined by 

= 2{D^e,.i,D^mitie,.i+D^G'l) for t S (t._i, t,], z = 1,2,..., n. 
It is easy to see that 


^ \Z^\^ dt < 4 max ^ 

l<2<'n .i—/ \ 




l|G^II/ds 


k \ 2=1 




Hence by the Davis inequality we have 


< E sup [ Z^ < 3E ( / V \Z^\^ ds] 
t^[0,T']Jo 0 ^ J 


1/2 


< ^E max ||DT'e,f-b36TVE||e*-i||Li-b36E / |lG.||f ds < oo. 
2 l<i<ra ^ Jo 

Using this, and estimates (13.311) and (13.431) . we obtain (hi) from (I3.44|l . 


□ 


By virtue of Remark 13.81 by using Sobolev’s theorem on embedding into Cb 
for m > d/2, we get the following corollary. 
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Corollary 3.11. (i) Let Assumptions 13.21 through 13.51 hold with m > d/2. Then 
for sufficiently small r there exists for all h and i? > 0 a unique solution = 

. 2; g HsuppCr} to ((3^ - (|3.10p . 

(ii) If Assumption [3i6] holds with some integer I > d/2 and Assumptions 13 .1 1 through 
13.51 hold with m > d/2 + 3, then 

(3.45) max E max \u^f^ — < Nt{1 +E/C^), 

0<i<n x^Gh 

with a constant N = N{K, C, H, d,T,l, Ag). 

(iii) If Assumption 13.61 holds with I > d/2 + 1 and Assumptions 13.ll through 157 
hold with m > (i/2 + 4, then 


Emax max < Nt{1 + E/C^) 

i<n x^Gh 


(3.46) 

with a constant N = N{K, C, H, d,T,l, Ag). 


Proof of Theorems \S. 5\ and \S.6] The solvability of (13.911 - (13.101) has already been 
discussed above, so we need only prove the estimates (13.121) and p.l5l) . A natural 
way to separate the errors of the different types of approximations would be to 
write 


i{x) - 


h,R,T 


(a;)| < \uriix) - u/:,{x )\-b {u/i^ix) - u'/:f{x)\ + {uffix) - 


h,Rf 


h,R,r 


{x)\. 


However, in such a decomposition we cannot directly estimate the middle term on 
the right-hand side, that is, the localization error of the finite difference equation. 
Therefore we introduce the classical solution of (CT)-(I^ with data 

® := = (Cr'01 Cr^i Crc, Cr(x, CrMi Cr/: Cr9)- 

Clearly the pair J), 2) satisfies Assumption 12.41 Also, as the finite difference co¬ 
efficients a, fa are compatible with the data (9 in the sense that Assumption 13.11 is 
satisfied, it follows that the coefficients fa-^, as defined in (I3.7ll - (I3.8L are com¬ 
patible with the data 9 in the same sense. Therefore in the decomposition 


{x) - •«'"(x)| < Mx) - uT{x)\ + \u/f{x) - Kf{x)\ + \u'/fix) - v^’^ix)] 


, 0,R 


0,R/ 


h,R( 


h,R/ 


h,R,T / 


the first term can be treated by Theorem 12.61 the second by Theorem 13.91 (c), and 
the third by Corollary [STT] (iii). Adding up the resulting errors we get the estimate 
(I3.12|) . and the dependence of the constants also follows from the invoked theorems. 
Similarly, for p.l5p we write 

\uTi{x) - Vi’^’"'{x)\ < \uriix) - (a;)| -b \u°f,^{x) - u'/:f{x)\ 


+ 


j =0 


|^./2^R(^)_^V2LR,.(^)|^ 


where and follow the same steps as above. 


Remark 3.12. As it can be easily seen from the last step of the proof. Assumption 
13.61 can be weakened to a-H61der continuity for any fixed a > 0, at the cost of 
lowering the rate from 1/2 to a A (1/2). 

To decrease the spatial regularity conditions, in particular, the term d/2, one 
would need the generalization of the results of [7], and subsequently, of Theorem 
13.101 to arbitrary Sobolev spaces W/f". Partial results in this direction can be found 
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